A Negative probability table problem of heating number in FENDL-3.1d ACE file

Konno, Chikara; Kwon, Saerom*

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.320 - 325, 2020/10

We found out that p-table data in the FENDL-3.1d ACE file included negative values for 33 nuclei. Thus, we studied why p-tables for heating number included negative. As a result, it was found out that partial KERMA factors became too large because the energy-balance was broken in the 33 nuclei and that FENDL-3.1d adopted kinematics KERMA factors. Then NJOY could not process adequately the 33 nuclei data, which led to negative p-tables for heating number. We prosed two solutions for this issue, produced new ACE files of FENDL-3.1d with the above two methods and confirmed that the new ACE files had no negative p-tables of the heating number.

Benchmarks of depletion and decay heat calculation between MENDEL and MARBLE

Yokoyama, Kenji; Lahaye, S.*

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.109 - 116, 2020/10

CEA/DEN/DM2S/SERMA and JAEA/NSEC are working on benchmarks for burnup, isotopic concentrations and decay heat calculations in the collaboration framework between both organisms. Both actors of this benchmark are independently developing their own simulation code systems for computing quantities of interest in nuclear fuel cycle domain: MENDEL in CEA and MARBLE in JAEA. The purpose of the benchmark is to verify each system by comparing both calculation results on specific applications. MENDEL uses a several solvers for the resolution of Bateman equation. Runge-Kutta method or Chebyshev Rational Approximation method (CRAM) are used for irradiation computations. An analytical solver can also be used for decay calculations. MARBLE can use Krylov subspace method or CRAM method. As the first phase of the benchmark, we compared the calculated results of decay heat and isotropic concentrations following by a Pu-239 fast fission pulse. We applied nuclear data from three libraries: (1) JEFF-3.1.1, (2) JENDL/DDF-2015 + JENDL/FPY-2011, and (3) ENDF/B-VII.1. Nuclear data and burnup chain were generated from these libraries independently on each system. We confirmed that the results for both systems were in very good agreement with each other. Numerical results were also compared to experimental data. As the second phase of the benchmark, we are proceeding with a burnup calculation benchmark of MENDEL and MARBLE using the nuclear data and burnup chain provided by ORLIBJ33, which is a set of cross-section data based on JENDL-3.3 for ORIGEN-2 code system. We will also compare with calculation results by the ORIGEN-2 code with ORLIBJ33. Since the series of ORLIB, that is, ORLIBJ32, ORLIBJ33, and ORLIBJ40, have been widely used especially in Japan for many years, the comparison with ORLIB is effective for confirming the performance of MENDEL and MARBLE.

An Investigation on the control rod homogenization method for next-generation fast reactor cores

Takino, Kazuo; Sugino, Kazuteru; Oki, Shigeo

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.92 - 96, 2020/10

Interactive in-situ steering and visualization of GPU-accelerated simulations using particle-based volume rendering

Kawamura, Takuma; Hasegawa, Yuta; Idomura, Yasuhiro

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.187 - 192, 2020/10

In order to realize the atmospheric dispersion prediction of pollutants, a fluid simulation by adaptive mesh refinement (AMR) optimized for GPU supercomputer has been developed, and interactive visualization and parameter steering of the simulation results are needed. In this study, we extend particle-based in-situ visualization method for structured grids into AMR, and enables in-situ steering of the simulation parameters by utilizing an in-situ control mechanism via files. By combining the developed method with plume dispersion simulation in urban areas running on a GPU platform, it was shown that human-in-the-loop pollution source search is possible without enormous parameter scanning.

Validation study of finite element thermal-hydraulics analysis code SPIRAL to a large-scale wire-wrapped fuel assembly at low flow rate condition

Yoshikawa, Ryuji; Imai, Yasutomo*; Kikuchi, Norihiro; Tanaka, Masaaki; Gerschenfeld, A.*

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.73 - 80, 2020/10

A finite element thermal-hydraulics simulation code SPIRAL has been developed in Japan Atomic Energy Agency (JAEA) to analyze the detailed thermal-hydraulics phenomena in a fuel assembly (FA) of Sodium-cooled Fast Reactors (SFRs). The numerical simulation of a large-scale sodium test for 91-pin bundle (GR91) at low flow rate condition was performed for the validation of SPIRAL with the hybrid k-e turbulence model to take into account the low Re number effect near the wall in the flow and temperature fields. Through the numerical simulation, specific velocity distribution affected by the buoyancy force was shown on the top of the heated region and the temperature distribution predicted by SPIRAL agreed with that measured in the experiment and the applicability of the SPIRAL to thermal-hydraulic evaluation of large-scale fuel assembly at low flow rate condition was indicated.

Image-based view selection for shape comparison of mode water regions in virtual reality spaces

Yano, Midori; Ito, Takayuki*; Tanaka, Yusuke*; Matsuoka, Daisuke*; Araki, Fumiaki*; Czauderna, T.*; Stephens, K.*

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.201 - 208, 2020/10

Visual analysis plays an important role in understanding and evaluating the climate models, their variables, and their outputs because complex processes are required for the tuning of the climate models. Virtual Reality (VR) technologies are effective for 3D visualization and have been recently employed for visual analysis of more various scientific data. However, researchers may have some problematic situations while using VR space depending on user operations and target geometries. Here, we propose an image-based view selection method to solve these situations and understand differences in ocean states between simulations and observations based on shapes of mode water regions. This view selection takes evaluation criteria for shape comparison of mode water regions generated from simulations and observations into account. This paper introduces two example cases applying this view selection and discusses on a degree of shape matching of mode water regions of each case.

Intuitive interactions for immersive data exploration of numerical simulation results

Tabata, Ginga*; Sakamoto, Naohisa*; Kawamura, Takuma

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.193 - 200, 2020/10

In interactive visualization/analysis, interactive operation of visualization results on the screen is an important technology that is directly linked to knowledge extraction. On the conventional 2D display, interactive operations such as moving the viewpoint and adjusting the visualization parameters were performed using the mouse. Then, the virtual reality (VR) technology such as CAVE system and head mounted display (HMD) is used for visualization and analysis, and the development of technology for interactively operating visualization result data in VR space has become one of the important issues in the visualization field. In this paper, we propose a technique to capture a real hand gesture by Leap Motion and intuitively change the viewpoint movement and visualization parameters in VR space by HMD. We asked four people to test the developed technology, and obtained a result that it was easy to change the viewpoint from the questionnaire on the operation feeling.

Ensemble wind simulations using a mesh-refined lattice Boltzmann method on GPU-accelerated systems

Hasegawa, Yuta; Onodera, Naoyuki; Idomura, Yasuhiro

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.236 - 242, 2020/10

The wind condition and the plume dispersion in urban areas are strongly affected by buildings and plants, which are hardly described in the conventional mesoscale simulations. To resolve this issue, we developed a GPU-based CFD code using a mesh-refined lattice Boltzmann method (LBM), which enables real-time plume dispersion simulations with a resolution of several meters. However, such high resolution simulations are highly turbulent and the time histories of the results are sensitive to various simulations conditions. In order to improve the reliability of such chaotic simulations, we developed an ensemble simulation approach, which enables a statistical estimation of the uncertainty. We examined the developed code against the field experiment JU2003 in Oklahoma City. In the comparison, the wind conditions showed good agreements, and the average values of the tracer gas concentration satisfied the factor 2 agreements between the ensemble simulation data and the experiment.

Molecular dynamics study of phosphorus migration in 5 grain boundary of -iron

Ebihara, Kenichi; Suzudo, Tomoaki

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.65 - 69, 2020/10

Phosphorus (P) is known as one of the elements which cause the grain boundary (GB) embrittlement in steels and its GB segregation is promoted by the increase of vacancies and self-interstitial atoms due to irradiation. Thus we have been developing the rate-theory model for estimating GB P segregation under several temperatures and irradiation conditions. Because the model does not include the trapping and de-trapping processes properly, however, the model cannot calculate GB P coverage which is measured by experiments. As for the de-trapping process, so far, we have considered the migration of a P atom in the GB region of 3 symmetrical tilt GB using molecular dynamics (MD). In the current study, we also simulated the P migration in 5 GB using MD and compared the result with that of 3. As a result, at 800K, it was found that a P atom cannot migrate in 5 without vacancies while a P atom can migrate between iron atoms in 3.

Machine-learning molecular dynamics study of thermal properties of CaF

Nakamura, Hiroki; Machida, Masahiko

Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.104 - 108, 2020/10